OPEN 3 ACCESS Freely available online 



•0-PLOS I o-^E 



Quantitative Estimation of Renal Function with fjfi 

Dynamic Contrast-Enhanced MRI Using a Modified cros^rk 
Two-Compartment Model 

Bin Chen\ Yudong Zhang^, Xiaojian Song^'^, Xiaoying Wang^'^*, Jue Zhang^'^% Jing Fang^'^ 

1 Academy for Advanced Interdisciplinary Studies, Pel<lng University, Beijing, China, 2 Department of Radiology, Peking University First Hospital, Beijing, China, 
3 Department of Electrical Engineering, Chengdu University of Information Technology, Chengdu, Sichuan, China, 4College of Engineering, Peking University, Beijing, 
China 



Abstract 

Objective: To establish a simple two-compartment model for glomerular filtration rate (GFR) and renal plasma flow (RPF) 
estimations by dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI). 

Materials and Mettiods: A total of eight New Zealand white rabbits were included in DCE-IVIRI. The two-compartment 
model was modified with the impulse residue function in this study. First, the reliability of GFR measurement of the 
proposed model was compared with other published models in IVlonte Carlo simulation at different noise levels. Then, 
functional parameters were estimated in six healthy rabbits to test the feasibility of the new model. Moreover, in order to 
investigate its validity of GFR estimation, two rabbits underwent acute ischemia surgical procedure in unilateral kidney 
before DCE-IVIRI, and pixel-wise measurements were implemented to detect the cortical GFR alterations between normal 
and abnormal kidneys. 

Results:lhe lowest variability of GFR and RPF measurements were found in the proposed model in the comparison. Mean 
GFR was 3.03±1.1 ml/min and mean RPF was 2.64±0.5 ml/g/min in normal animals, which were in good agreement with 
the published values. Moreover, large GFR decline was found in dysfunction kidneys comparing to the contralateral control 
group. 

Conclusion: Results in our study demonstrate that measurement of renal kinetic parameters based on the proposed model 
is feasible and it has the ability to discriminate GFR changes in healthy and diseased kidneys. 
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Introduction 

Glomerular filtration rate (GFR) is the most critical functional 
parameter of the kidney [1,2]. InuKn clearance is widely 
considered the gold standard in quantification of GFR. However, 
it is not routinely being used in the clinical setting, since it requires 
multiple blood and urine samples over a period of several hours 
and could not obtain single-kidney GFR. Traditional radioactive 
labeled methods are invasive and result in radiation burden. 
Therefore, it is important to establish a simple new method for 
quantitative measurements of GFR to monitor renal function. 

Recently, on the basis of dynamic contrast-enhanced magnetic 
resonance imaging (DCE-MRI), functional parameters could be 
estimated in a wide range of applications including breast imaging 
[3], brain imaging [4], and abdominal imaging [5]. In order to 
overcome some mentioned limitations of GFR measurement, 
several compartment models based on DCE-MRI have been 
proposed. Patlak-Rutland method was introduced to evaluate 
glomerular filtration rate, neglecting the outflow of contrast agent 
from the tubular compartment [6] . Under this assumption, it was 



often appUed to the whole kidney to ensure the tracer stays in the 
region of interest (ROI) [7,8], which, however, could not directly 
reflect GFR since the glomerular filtration predominandy exist in 
renal cortex [9-12]. 

Subsequentiy, a two-compartment model with constant disper- 
sion of contrast agent and the consideration of outflow was 
proposed [13]. However, lacking of accurate tubule parameters 
estimation was supported and underestimated GFR was found. 
Later, a similar model was employed to measure GFR in patients 
by using the whole kidney ROI, which was just simplified by 
ignoring the dispersion effect in the original two-compartment 
model [14]. Then, a separable compartment model, of which the 
time delay of contrast was assumed to be zero [15], provided 
estimated parameters which were analogues to those of the 
original two-compartment model [13]. 

Besides these models, three-compartment models were also 
proposed for GFR estimations [16,17]. Nevertheless, the variabil- 
ity of GFR measurement was more pronounced than that of two- 
compartment models [18]. In addition, the assessment of GFR was 
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sensitive to tlie segmentation of cortex and medulla, which may 
result in an unstable measurement [19]- The increase of 
compartment number would result in an expense of system 
complexity, and a simple cortical compartment model is alterna- 
tive because the glomerular filtration occurs in the cortex. We 
borrowed Zhang et al's idea of using impulse residue function to 
improve the robustness in GFR measurements [17], and combine 
it into a two-compartment model. 

In this study, a modified two-compartment model was proposed 
by introducing the impulse residue function to obtain effective 
estimations of GFR and renal plasma flow (RPF) from DCE-MRI. 
The advantage of the new model in GFR or RPF measurements 
over other published models was investigated in Monte Carlo 
simulation under different noise levels. Then, quantitative 
estimations were performed in healthy rabbits by using our new 
model. Furthermore, the proposed model was employed to 
measure GFR changes in rabbits with unilateral ischaemic acute 
kidney injury (AKJ) to test its validity. Pixel-wise calculation was 
performed and the cortical GFR results were compared with its 
contralateral kidney. The major focus of this study was put on the 
sufficient robustness in kinetic parameters estimation and the 
ability to discriminate healthy and diseased kidneys. 

Materials and Methods 

Animals 

This study was approved by the Ethics Committee on animal 
care and use of Peking University First Hospital (Ethic number: 
10/235). Eight male New Zealand white rabbits (weighing from 
2.9 kg to 3.5 kg) underwent this experiment (ii = 6 for normal 
kidney experiment; n = 2 for dysfunction kidney experiment). AU 
the rabbits were housed in individual cages at room temperature, 
and free fed with standard feed and tap water, and stopped 
12 hours before experiments. Before DCE-MRI, pentobarbital 
sodium (0.5 ml/kg body mass) was injected through the marginal 
ear vein with a 24-gauge catheter for anesthetization. All the 
animals were placed in a supine position in a fixed device to limit 
abdominal motion during scans. Heart rate was continuously 
monitored by using a photopulse sensor of the MR scanner during 
acquisitions. 

In the experiment of dysfunction kidneys, two rabbits (2.9 kg 
and 3.3 kg) were implemented a unilateral renal ligation surgery 
procedure after anesthetization. After a midline incision, the aorta 
of left kidney was clamped for 45 minutes to induce ischaemic 
acute kidney injury. The right kidney was normal and regarded as 
the control group. After surgery, we sew up the incision with 
sewing wires and give analgesic treatment to minimize suffering. 



The animals' body temperature was maintained about 38° by 
using a hotplate during surgery and exposed to an infrared light in 
the box before MR scans. Two days later, the DCE-MRI 
experiment was implemented. These two injured rabbits (acute 
ischemia in the left kidney) were scanned together with all the 
other six normal rabbits. 

MRI 

The abdomen dynamic images were acquired on a whole body 
3.0 T MR scanner (Signa Excite; GE Medical Systems, Milwaukee, 
WI, USA) with TORSOPA coil. Before each DCE-MRI scan, Tl 
measurement was conducted firstly with a three-dimensional spoiled 
gradient-recalled echo sequence. Tl mapping was measured by 
using variable fUp angles method [20] , and the imaging parameters 
are as follows: TR = 6.2 msec, TE = 2.9 msec, flip angle = [3,9,20] °, 
slices =16, slice thickness = 4.0 mm, matrix = 256x256, FOV = 
200 X 200 mm . Then, DCE-MRI scans were performed by a three- 
dimensional fast spoiled gradient recalled echo sequence with the 
following parameters: TR= 3.3 msec, TE= 1.3 msec, flip an- 
gle =15°, FOV= 160x160 mm^, ASSET = 2, acquisition matrix 
was 128x128 and interpolated to 256x256, slice thick- 
ness =4.0 mm, slice number = 16, bandwidth was 488 Hz/pixel, 
and acquisition time was 3.0 s/frame. Five precontrast frames were 
obtained before bolus injection, then, an administration of 
0.05 mmol/kg of Gd-DTPA (Omniscan; GE Healthcare Ireland, 
IDA Business Park, CarrigtohiU, Co.Cork) was performed with a 
venous cannula at a rate of 2 ml/ s. 5 ml of saline was immediately 
flushed in and totally 100 frames were obtained in about five 
minutes. 

Data Analysis 

Renal parenchyma and the aorta were automatically segmented 
from the surrounding tissue with a Level-Set framework [21]. The 
parameters used in segmentation were: alpha = 0.04, a parameter 
that controls the weight of smoothing item of the image features; 
and iteration was set to 60, such values were sufficient for detection 
of kidney outline in this experiment. Afterwards, dynamic images 
of the segmented regions were registrated to reduce motion [22]. 
Then, a slice which covering the largest possible parenchyma 
during corticomeduUary phase was used for cortical ROI drawing 
to obtain the tissue signal intensity curves. The size of ROIs for all 
the kidneys was 590 pixels in average. The arterial input function 
(AIF) was determined by drawing ROI in the slice that a branch of 
renal artery was obviously seen, that is, the ROI was placed within 
the aorta distal to the branch of the renal artery. To improve the 
robustness of AIF and reduce the inflow artifacts, one more slice 



A„(t) 



A 



R.(t) 



o,{') 



o 



R/t) 



-ignllRF, 



Figure 1. Schematic diagram of the modified two-compartment model with impulse residue function for glomerular filtration. A Is 

the vascular compartment Includes Intrarenal arteries and glomerular vessels, and T Is the tubules compartment. The retention function and Rj In 
compartment A and T (represented In solid arrow within the box) are the convolution of the Input and each Impulse residue function. Dashed lines 
denote the outflow of each compartment and the outflow of compartment A partially flows Into compartment T. 
dol:1 0.1 371/journal.pone.01 05087.g001 
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Figure 2. Concentration curves in Monte Carlo simulation for a.) the artificial AIF and b.) the tissue curve which is generated with 
the known initial parameters and added with 5% level noise. 

doi:1 0.1 371/journal.pone.01 05087.g002 



above was also selected for averaging. A 3 x 3 pixel size ROI was 
used. Furthermore, the tail of the AIF was fit to a biexponential 
decay to reduce respiratory motion related noise [23]. 

Then, signal intensity curves were converted into gadolinium 
concentrations with the Tl values from Tl mapping [23], under 
the assumption that a linear relationship between relaxation rate 
changes and [Gd] concentrations according to the following 
equation: 

l/Tlpost = l/Tlp„ + ri*[Gd] (1) 

where, 1/Tlp„. and l/Tlpo,t are the contrast relaxation rates 
before and after bolus injection, ri is the specific Tl relaxivity at 
physiological temperatures in plasma (4.1 L s mmole ) 
[24,25]. 

The New Compartment Model 

The proposed model in this study describes two compartments: 
the intrarenal arteries and glomerular vessels (A) and the renal 
tubules in cortex (T), with tracer flowing from compartment A into 



compartment T, shown in Fig. 1 . The concentration in descending 
aorta is expressed as the arterial input Ao(t), and the concentration 
in cortex, C(.ortcx(t) is contributed by the retention in compartment 
A and T over time. The arterial input A()(t) is converted into 
plasma concentration Ap(t) by dividing by (1-Hct), where 
Hct = 0.45 is the hematocrit for rabbits. 

While contrast passes through the kidney, an administration of 
impulse residue function (IRF) can be useful to describe the 
characteristics of contrast agent distribution in each compartment 
[26]. In our study, considering the dispersion effect and transit 
time of contrast, piecewise-exponential impulse residue function is 
introduced to determine its response to the idealized bolus 
injection. For an ideal instantaneous unit bolus injection, it is 
actually a time-enhancement function [27]. The contrast is 
considered to be a bolus injection into the kidney and equals 1 
(t<T;), then an exponential fall with constant rate m; starts when 
the contrast is washing out (t>T;), and the parameter T; is the 
minimal transit time: 
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Figure 3. Monte Carlo simulation results for each compartment model in estimating a.) GFR and b.) RPF, only three models could 
extract RPF. Lowest variability of estimated GFR and RPF (red dotted line) are found in the proposed two-compartment model at 2%, 3%, 5%, 10% 
and 15% noise, respectively. 
doi:1 0.1 371/journal.pone.01 05087.g003 
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Figure 4. A representative image from DCE-MRI acquisition, a.) showing the region of interest after segmentation; b.) the manually drawn 
cortical ROI (white line); c.) the image in which a branch of renal artery is clearly seen is used for aortic ROI drawing; and d.) the concentration curves 
of each corresponding ROI. The original aortic input function (black line) fits its tail with a biexponential method (shown in red line). The blue line is 
the cortical concentration curve. 
doi:10.1371/journal.pone.0105087.g004 



imit)-- 



1 ?<T, 



(2) 



where the IRFi(t) is the impulse residue function and i represents 
A and T. 

Witli the proposed model (Fig. 1), the residue amount of 
contrast in each compartment after a bolus injection is presented 
as such a convolution of the input and its corresponding impulse 
residue function: 



CA(t) = RPF-Ap{t)®IRFA{t) 
Crit) = GFR- OA{t)®IRFT(t) 



(3) 



where ® express the convolution, RPF is the renal plasma flow, 
GFR is the glomerular filtration rate, and Ci(t) is the contrast 
concentration in each compartment. 



In renal system, we defined the retention function Ra and Rx to 
illustrate the changes of contrast for each compartment after the 
convolution of IRF;. We assumed the unit input for compartment 
A is 1, then the retention function of A equals IRFA(t), however, 
the unit input changes in compartment T. Because the input for T 
is from the output of compartment A, the retention function for T 
is the convolution of OA(t) and IRFr(t), then equation 3 becomes: 



where 



CA(t) = RPF-A„{t)®RA 
CT{t) = GFR-Ap(t)®RT 



RA=IRFA(t) 

1 t<-lA 



(4) 




Time (s) 

Figure 5. Typical retention function curves for vascular 
compartment (Ra, red), tubule compartment (Rt, green) and 
the cortex region (Rcortex, blue). 

doi:1 0.1 371/journal.pone.Ol 05087.g005 



Rr = m- ^''^]®IRFT{t) 

r 0 t<TA 

= }[\-e-'"A<~'-^A)] TA<t<TA+TT 
[[El-e-'"A('-U)] t>TA+TT 

3{t) is the unit input and then cortical retention function written 

as: Rcortex =Ra+Rt- 

The mean transit time (MTT) was defined as the area under the 
retention curves according to the equation: MTT-, = Rj(j)dt. 

Thus, the amount of contrast remained in cortex region is given 
by: 
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Figure 6. GFR mapping of a.) tKie normal kidney (right kidney) and b.) tKie acute iscKiemia kidney (left kidney) by using the new 
model. The left kidney of the rabbit is under surgical ligation for totally 45 minutes, and lower GFR values are clearly observed in the cortex and 
outer stripes of the outer medulla regions while corresponding regions are high in the control group. 
doi:10.1371/journal.pone.0105087.g006 



\ (5) 



\-Hct 



{RPF-GFR)e-'"A(i-'A> + GFREi t>XA + tT 



where Ei = -t „,^_^^ , C„,rt,,x(t) is the 

concentration in cortex. Parameters were fitted by using the 
nonlinear least-squares algorithm. 

Monte Carlo Simulation 

To evaluate the reliability of the new model compared with 
other published models in estimating GFR in rabbit kidney, 
Monte Carlo simulations were conducted. Artificial aortic input 
function used in this simulation was obtained by smoothing an AIF 
curve with temporal resolution of 3 second. Smoothing was carried 
out by a gamma variante function [28]. 

Totally six compartmental models were included with the same 
preset GFR values in this simulation: Patlak-Rutland method 
(GFR = 3.0 ml/min, v„ = 2.4), PR [6]; a two- compartment model 
(GFR=3.0 ml/min, k„„, = 0.08, 4 = 0.2, d = 0.42), 2CD [13]; a 
simplified two-compartment model (GFR = 3.0 ml/min, k„^, = 
0.08, 4 = 0.2), 2C [14]; a separable compartment model (GFR = 
3.0 ml/min, T,= l/k„„t =12.5,V,,= f, =0.2,Tp= d =0.42), SP 
[15]; a three-compartment model (GFR = 3.0 ml/min, RPF = 
3.0ml/g/min, fp = 0.18, w,j = 0.2, Wa_na = 0.08, Wp = 0.18), 3C 
[16]; and the proposed two compartment model with impulse 
residue function in this study(GFR = 3.0 ml/min, RPF = 3.0 ml/ 
g/min, mA = 0.2, t:a= 3.2 s, mx = 0.1) and Tt was assumed to be 
0, 2C-IRF. Then, the tissue curves were generated from their 
corresponding initial parameters and the artificial AIF. Consider- 
ing the impact of noise in DCE-MRI scans, different levels of noise 
were reintroduced to each concentration curve in Monte Carlo 
simulation. Guassian noise with zero mean and SD equals to 2%, 
3%, 5%, 10%, and 15% of the mean magnitude of each 
gadolinium concentration curve were used. These noise were 
generated by using randn(l,N) function in Matiab and be added to 
construct noisy tissue curves. AH these models were used to fit its 
corresponding noisy data for which true values of these parameters 
were known. 2000 Monte Carlo trials were conducted, and the 
variability of each fitted parameter was obtained. In order to assess 



the variability of the estimated parameter, the coefficient of 
variation (CV) was used according to the equation: CV = SD/ 
mean. The bias in estimated parameters can be given by the 
difference between mean value of 2000 simulations and the actual 
values. 

In-vivo Experiment 

After testing the reUabUity and bias of the proposed model in 
parameters estimation and comparing the precision of GFR 
measurement with other published models, the extracted concen- 
tration curves from DCE-MRI data were used for the measure- 
ments of kinetic parameters with 2C-IRF model in normal rabbit 
kidneys. Nonlinear least squares fitting was implemented with all 
the data sets by using Levenberg-Marquardt algorithm. To assess 
the goodness of fit, coefficient of determination, denoted R , was 
also calculated. 

Acute kidney injury causes significant renal dysfunction in terms 
of decreased glomerular filtration rate [29] . In order to investigate 
the validity of the proposed model to detect GFR alterations in 
such an abnormal condition, pixel-wise calculation of GFR were 
measured with the dysfunction kidney, as well as the contralateral 
kidney. To further quantify the cortical mean GFR in the 
abnormal kidneys, two-certified radiologists with at least 5 years of 
experience in renal DCE MR imaging, who were blinded to the 
experiment, draw the cortical ROI on the Tl -weighted images 
and the mean GFR values from its corresponding mapping image 
were obtained. 

Histology 

Kidneys in dysfunction experiment were fixed in 10% neutral 
buffered formalin and embedded in paraffin for light microscopic 
study. Kidneys were sectioned into S-^im slides and stained for 
histology with hematoxylin-eosin. One experienced pathologist, 
who was blind to which experimental group the samples belonged, 
reviewed histological findings. 

Statistical Analysis 

Results of the quantitative measurements are expressed as mean 
+ SD. Paired t-test was performed and a P value<0.05 was 
considered statistically significant. Coefficient of determination 
(R^) is calculated to assess the goodness of fit. All analysis are 
implemented in Matiab (MathWorks, Natick, MA). 
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Figure 7. Histology sections of a.) the normal kidney (H&E, x100) and the acute ischemia kidney with histological findings of b.) 
tubular dilatation (x100), c.) cast deposition (open arrow, x100) and d.) cell necrosis (black arrow, x200). 

doi:10.1371/joumal.pone.0105087.g007 



Table 2. Pixel-wise estlnnatlon of GFR for dysfunction kidneys by using the proposed two-compartment model. 







Cortical 


Region 


Kidneys 


pixel GFR 


mean GFR 




(ml/min) 


(ml/min) 


Injury 1. 


0.0037 


2.27' 


Control 1. 


0.0080 


4.55 


Injury 2. 


0.0042 


2.78' 


Control 2. 


0.0068 


4.73 



'pronounced reduction was found compared to the control group. 
doi:l 0.1 371 /journal.pone.Ol 05087.t002 
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Table 3. Estimated results of GFR for each two-compartnnental nnodel. 





n = 6 


2C-IRF 


2CD 


2C 


SP 


PR 


GFR{ml/min) 


3.03 


1.29 


1.55 


1.80 


0.84 


(±SD) 


1.1 


0.7 


1.1 


0.6 


0.5 



doi:l 0.1 371 /journal.pone.Ol 05087.t003 



Results 

Monte Carlo Simulation 

In Monte Carlo simulation, the artificial AIF is shown in 
Fig. 2a, and 2%, 3%, 5%, 10% and 15% noise were reintroduced 
into the concentration curves at all the time points, shown in 
Fig. 2b with 5% noise for instance. 

Compared with other five published models, the smallest CV 
values of estimated GFR were found in 2C-IRF model, seen in 
Fig. 3a. Higher variability was found in the other compartment 
models in GFR estimate and pronounced impact existed in 3C 
model when noise level was above 5%. The bias of estimated GFR 
was 0.01%-4.0% for 2C-IRF, 0.2-11.3% for 2C, 2.6%-9.8% for 
2CD, 8.4%-14.7% for SP, 1.1%-17.1% for 3C, and 0.2%-1.7% 
for PR. For RPF Monte Carlo simulation, smallest CV values were 
also found in 2C-IRF model, ranging from 0.2% to 2.0%, seen in 
Fig. 3b. The other models have a bigger variability of more than 
10% when the noise level is above 5%. The bias of estimated RPF 
of 2CTRF model was 0.05% to 0.3%, much lower than those of 
other models. 

The estimated GFR and RPF results are robust under noise 
conditions by using 2C-IRF model. Small variability of estimated 
RPF was found with CV values ranging from 0.2% to 2.0%, while 
CV values of GFR were 4.2% to 10.8% under 10% noise. The 
bias of the estimated GFR, RPF from the new model was low with 
a range of 0.0%-0.6% for different noise levels, excluding one 
larger bias of 4.0% for GFR at 15% noise level. 

In-vivo Experiment 

In this study, all the rabbits were analyzed with the 2C-IRF 
model for kinetic parameter estimations. A segmentation result is 
shown in Fig. 4a. After registration, cortical ROI was manually 
drawn (shown in Fig. 4b) and slice that with branch of renal artery 
(white arrow in Fig. 4c) was selected for aortic ROI drawing. 

Table 4. RPF values in literature studies. 



Authors 


year 


ml/g/min 


method 


Ott CE et al. [31] 


1979 


3.2±0.3 {n = 13) 


inulin clearance 


Hermoye L et a!.*" [30] 


2004 


1.3±0.4 (n = 6) 


deconvoiution 


Winter JD et al. [22] 


2011 


3.28±0.59 {n = 5) 


ASf' based 






2.98±0.60 (n = 6) 


DCE*" based 


Zhang Y et al. [32] 


2012 


3.20±0.67 (n = 9) 


ASL based 


Zimmer F et al.*'' [33] 


2013 


4.1 6± 1.24 (n = 6) 


ASL based 






5.42±0.85 (n = 6) 


DCE based 



Unit Is ml/g/min. 

*^unit converted to ml/g/min with 1 g/ml density. 

*'^rat models were used in RPF measurements. Five of the six rats had a unilateral ischaemic AKI and total seven healthy kidneys were included in the calculation of 
mean RPF. 

•'arterial spin labeling (ASL) MRI. 
**^dynamic contrast-enhanced (DCE) MRI. 
doi:l 0.1 371/journal.pone.Ol 05087.1004 



Representative concentration curves derived from corresponding 
ROIs are shown in Fig. 4d. 

Estimation results of all the kidneys are listed in Table 1 . GFR 
obtained from the new model was 3.03±1.1 ml/min, RPF was 
2.64±0.5 ml/g/min, vascular mean transit time MTT^, was 
5.6±0.6 s, tubule mean transit time MTT^ was 15.1 ±8.8 s and 
the mean transit time of the kidney MTTk was 20.7±8.7 s in 
average. No significant differences were found for GFR, RPF, 
MTTa, MTTx or MTTk when comparing the left and right 
kidney groups (p>0.05 for all). Satisfactory goodness of fit was 
found for all the cases, with R^ ranged from 0.81 to 0.97. Typical 
retention curves are shown in Fig. 5. 

In dysfunction kidney experiments, large decrease of GFR was 
found in the cortical region of the left kidneys, shown in Fig. 6b, 
while the right kidneys (control group. Fig. 6a) remain normal. 
Pixel-wise calculation was implemented and the GFR mapping 
shows the differences between the dysfunction kidney and the 
contralateral normal kidney, which are in good concordance with 
the cortical mean values that listed in Table 2. The changes of the 
kidney with acute injury are particularly visible on histology 
compared to the normal kidney (Fig. 7a). Histological findings 
show main tubular dilatation (Fig. 7b), cast deposition (open 
arrow. Fig. 7c) and cell necrosis (black arrow, Fig.7d) in the 
injured kidney. 

Discussion 

In present study, a modified two-compartment model with 
impulse residue function is implemented, and Monte Carlo 
simulation results indicate the reliability of the proposed model 
in parameters estimates. The main result of this study is that the 
GFR and RPF measurements are in close agreement with 
literature results in rabbits. Moreover, the new model is valid in 
detecting GFR alterations in diseased kidneys. The GFR value of 
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normal kidneys in our study is 3.03 ±1.1 ml/ min (n = 6), which is 
comparable to a mean value of 2.9±1.0 ml/min measured in 
rabbits using a deconvolution method [30] and model derived 
GFR of 2.3±1.0 ml/min with original two-compartment model 
and 2. 8± 1.2 ml/min with Padak-Rudand method reported in 
rabbits study (n= 10) [13]. Here, GFR results of each compart- 
ment model for in-vivo experiments are presented in Table 3. The 
close agreement with literature results for the proposed model is 
found. 

Moreover, the RPF measured in this study is 2.64±0.5 ml/g/ 
min, which is in close agreement with literature reports [22,30- 
33], shown in Table 4. Our results resemble the renal perfusion 
reported in rabbits and larger than the deconvolution method 
(Hermoye L et al.). The perfusion results in rats (Zimmer F et al., 
4.16±1.24 ml/g/min based on arterial spin labeling method and 
5.42±0.85 ml/g/min based on DCE method) are larger than the 
results in our study (2.64±0.,5 ml/g/min), this discrepancy 
underscores physiological differences across species. The close 
agreement with literature and the vaHd assessments of discrimi- 
nating healthy and dysfunction kidneys indicate that the proposed 
2C-IRF model provides a useful method for quantitative 
measurements in kidneys, and the dilference in accuracy is 
relevant. Our results show that the new model is feasible and not 
only achieves reliable measurements of renal function, but also be 
capable in detecting of GFR alterations in dysfunction kidneys. 

Using the predefined impulse residue function, which considers 
the transit time and dispersion effect of contrast agent in kidney, 
our new model could characterize the distribution of contrast and 
elucidate physiological meaning better, as well as improves the 
model reliability with reduced bias in measurements of parame- 
ters. A similar use of impulse residue function approach has been 
reported [17]. However, in order to increase the accuracy of the 
estimated parameters, they introduced more compartments, which 
results in the expense of system complexity. 

In this study, the smallest variability of estimated GFR and RPF 
is found in 2C-IRF model, which illustrates that it is not sensitive 
to noise and sufficiently roljust in GFR or RPF measurements 
from DCE-MRI. Usually, deconvolution method is clinically used 

References 

1. Hayashi T, Nitta K, Hatano M, Nakauchi M, Nihci H (1999) The serum 
cystatin C concentration measured by particle-enhanced immunoncphelomctry 
is well correlated with inulin clearance in patients with various types of 
glomerulonephritis. Nephron 82:90-92. 

2. Craig AJ, Britten A, Heenan SD, Irwin AG (201 1) Signiiicant differences when 
using MDRD for GFR estimation compared to radionuclide measured 
clearance. Eur Radiol 21:2211-2217. 

3. FlucldgerJU, Schabel MC, DibeUa EV (2012) The effect of temporal sampling 
on quantitative pharmacoldnetic and three-time-point analysis of breast DCE- 
MRI. Magn Reson Imaging .^80:934-943. 

4. Pauliah M, Saxena V, Maris M, Husain N, Rathore RK, et al. (2007) Improved 
T(l)-weighted dynamic contrast-enhanced MRI to probe mierovascularity and 
heterogeneity of human glioma. Magn Reson Imaging 2.'i: 1292— 1299. 

5. Khalffa F, Abou El-Ghar M, AbdoUahi B, Frieboes HB, El-Diasty T, et al. 
(2013) A comprehensive non-invasive framework for automated evaluation of 
acute renal transplant rejection using DCE-MRI. NMR Biomed 26:1460-1470. 
doi:10.1002/nbm.2977. 

6. Hackstein N, Heckrodt J, Rau WS (2003) Measurement of single-kidney 
glomerular filtration rate using a contrast-enhanced dynamic gradient-echo 
sequence and the Rutland-Patlak plot technique. J Magn Reson Imaging 
18:714-72."). 

7. Hackstein N, Bauer J, Hauck EVV, Ludwig M, Kramer HJ, et al. (2003) 
Measuring single-kidnev glomerular filtration rate on single-detector helical CT 
using a Iwo-poinI Patlak plot technique in patients with increased interstitial 
space. AJR,\mJ Roentgenol 181:147-156. 

8. Tsushima Y, Blomley MJ, Kusano S, Endo K (1999) Use of contrast-enhanced 
computed tomography to measure clearance per unit renal volume: a novel 
measurement of renal function and fractional vascular volume. Am J Kidney 
Dis 33:754-760. 

9. Woff GL (1999) Using enhanced computed tomography to measure renal 
function and fractional vascular volume. Am J Kidney Dis 33:804—806. 



to determine renal transit time, however, it continues to exhibit 
noise in the deconvolved curves. By using the impulse residue 
function in our model, the vascular mean transit time is 5.6±0.6 s, 
which is similar to the result (6.8±1.7 s) in a previous study [30], 
and we also generated the tubule mean transit time (15.1 ±8.8 s) 
and kidney mean transit time (20.7±8.7 s). 

In DCE-MRI scans, administration of gadolinium contrast is 
usually used to enhance the signal. Gd contrast agents are rapidly 
cleared with a half-life of about 2 h in normal kidneys, however, it 
would exceed in patient with dysfunction kidneys. Thus, the 
retained and subsequent retention of Gd contrast would activate 
illness known as the nephrogenic systemic fibrosis (NSF) disease 
[34,35]. Researchers demonstrated Gd contrast possibly plays the 
triggering role in the development of NSF [36]. Thus, higher dose 
of Gd contrast would more easily lead to NSF disease. In our 
study, low dose of contrast about 0.05 mmol/kg was used and this 
may reduce the risk of NSF. 

There are several limitations in our study. First, the population 
of the dysfunction kidneys is small. Second, the arterial input 
function may be affected by the inflow eflects in the aorta. Last, 
the validity of the new model should be tested in human kidneys 
before clinical utility. 

In conclusion, our new model with the introduction of impulse 
residue function is feasible and suitable to estimate important renal 
functional parameters, and has the ability to discriminate GFR 
changes in healthy and diseased kidneys. 

Supporting Information 

Checklist SI ARRIVE Checklist. 

(DOC) 

Author Contributions 

Conceived and designed the experiments: BC YZ JZ. Performed the 
experiments: BC YZ XS XW. Analyzed the data: BC YZ XS JZ JF. 
Contributed reageiits/materials/analysis tools: BC JZ XWJF. Wrote the 
paper: BC. Interpretation of data; XS XWJF. Revising the article: BG YZ 
XS XWJZJF. Final approval: BC YZ XS XWJZ JF. 



10. Krier JD, Ritman EL, Bajzer Z, Romero JC, Lerman A, et al. (2001) 
Noninvasive measurement of concurrent single-kidney perfusion, glomerular 
filtration, and tubular function. Am J Physiol Renal Physiol 281: F630-F638. 

11. Laurent D, Poiricr K, Wasvary J, Rudin M (2002) Effect of essential 
hypertension on kidney lunetion as measured in rat by dynamic MRI. Magn 
Reson Med 47:127-134. 

12. Miles KA, Lcggctt DA, Bennett GA (1999) CT derived Patiak images of die 
human kidney. Br J Radiol 72:153-158. 

13. Annet L, Hermoye L, Peeters F, Jamar F, DehouxJP, et al. (2004) Glomerular 
filtration rate: assessment with dynamic contrast-enhanced MRI and a cortical- 
compartment model in the rabbit kidney. J Magn Reson Imaging 20:843-849. 

14. Buckley DL, Shurrab AE, Cheung CM, Jones AP, Mamtora H, et al. (2006) 
Measurement of single kidney function using dynamic contrast-enhanced MRI: 
comparison of two models in human subjects. J Magn Reson Imaging 24: 1117— 
1123. 

15. Sourbron SP, Miehacly HJ, Reiser MF, Schoenberg SO (2008) MRI- 
measiu-ement of perfusion and glomerular filtration in the human kidney with 
a separable compartment model. Invest Radiol 43:40-48. 

16. Lee VS, Rusinck H, Bokachcva L, Huang AJ, Oesingmann N, et al. (2007) 
Renal function measurements Irom MR renography and a simplified multi- 
compartmcntal model. Am J Physiol Renal Physiol 292:F1548-F1559. 

17. Zhang JL, Rusinck H, Bokachcva L, Lerman LO, Chen Q, et al. (2008) 
Functional assessment of the kidney from magnetic resonance and computed 
tomography renography: impulse retention approach to a multicompartment 
modd. Magn Reson Med 59:278-288. 

18. Bokacheva L, Rusinek H, Zhang JL, Chen Q, Lee VS (2009) Estimates of 
glomerular filtration rate from MR renography and tracer kinetic models. 
J Magn Reson Imaging 29:371-382. 

19. Rusinek H, Boykov Y, Kaur M, Wong S, Bokachcva L, et al. (2007) 
Performance of an automated segmentation algorithm for 3D MR renography. 
Magn Reson Med 57:1159-1167. 



PLOS ONE I www.plosone.org 



9 



August 2014 I Volume 9 | Issue 8 | e105087 



Quantitative Renal Function Measurement 



20. Cheng HL, Wright GA (2006) Rapid high-resolution T(l) mapping by variable 
flip angles: accurate and precise measurements in the presenceof radiofrequency 
field inhomogeneity. Magn Reson Med 55:566-574. 

21. Vese LA, Chan TF (2002) A multiphase level set fi^amework for image 
segmentation using the Mumford and Shah model. Int J Comput Vis 50:271— 
293. 

22. Lv D, Zhuang J, Chen H, Wang J, Xu Y, et al. (2008) Dynamic Contrast- 
Enhanced Magnetic Resonance Images of the Kidney. IEEE Eng Med Biol Mag 

27:36-41. 

23. Winter JD, St Lawrence KS, Cheng HL (2011) Qiiantification of renal 
perfusion: comparison ot arterial spin labeling and dynamic contrast-enhanced 
MRI. J Magn Reson Imaging 34:608-615. 

24. Peeters F, Annet L, Hermoye L, Van Beers BE (2004) Inflow correction of 
hepatic perfusion measurements using Tl -weighted, fast gradient-echo, contrast- 
enhanced MRI. Magn Reson Med 51:710-717. 

25. Shuter B, Wang SC, Roche J, Briggs G, PopeJM (1998) Relaxivity of Gd-EOB- 
DTPA in the normal and biliary obstructed guinea pig. J Magn Reson Imaging 
8:853-861. 

26. Miles KA, Griffiths MR (2003) Perfusion CT: a worthwhile enhancement? 
Br J Radiol 76:220-231. 

27. Lawson RS (1999) Application of mathematical methods in dynamic nuclear 
medicine studies. Phys Med Biol 44:57-98. 

28. de Priester JA, den Boer JA, Christiaans MH, Kessels AG, Giele EL, et al. (2003) 
Automated quantitative evaluation of diseased and nondiseased renal transplants 
with MR renography. J Magn Reson Imaging 17:95—103. 



29. Joyce M, Kelly C, Winter D, Chen G, Leahy A, et al. (2001) Pravastatin, a 3- 
hydroxy-3-methylglutaryl coenzyme A reductase inhibitor, attenuates renal 
injury in an experimental model of ischcmia-reperfusion. J Surg Res 101:79—84. 

30. Hermoye L, Annet L, Lemmerling P, Peeters F, Jamar E, et al. (2004) 
Calculation of the renal perfusion and glomerular filtration rate from the renal 
impulse response obtained with MRI. Magn Reson Med 51:1017-1025. 

31. Ott CE, Vaxi RC (1979) Renal autoregulation of blood flow and filtration rate in 
die rabbit. Am J Physiol 237:F479-F482. 

32. Zhang Y, Wang J, Yang X, Wang X, Zhang J, et al. (2012) The serial effect of 
iodinated contrast media on renal hemodynamics and oxygenation as evaluated 
by ASL and BOLD MRI. Contrast Media Mol Imaging 7:418-425. 

33. Zimmer E, Zollner E(}, Hoeger S, Klotz S, Tsagogiorgas CJ, et al. (2013) 
Quantitative renal perfusion measurements in a rat model of acute kidney injury 
at 3T: testing inter- and intramethodical significance iil' ASL £md DCE-MRI. 
PLoS One 8:c53849. doi: 10.1371/journal.ponc.0053849. 

34. Sadowski EA, Bennett LK, Chan MR, Wentland AL, Garrett AL, et al. (2007) 
Nephrogenic systemic fibrosis: risk factors and incidence estimation. Radiology 
243:148-157. 

35. Kuo PH (2008) NSF-active and NSF-inert species of gadolinium: mechanistic 
and clinical implications- AJR Am J Roentgenol 191:1861—1863. 

36. Grobner T (2006) Gadolinium- -a specific trigger for the development of 
nephrogenic fibrosing dermopathy and nephrogenic systemic fibrosis? Nephrol 
Dial Transplant 21:1104-1108. 



PLOS ONE I www.plosone.org 



10 



August 2014 I Volume 9 | Issue 8 | e105087 



